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Introduction 

Wear protective coatings have been developed to protect the contacting surfaces from 
friction and wear. The design architecture and composition of the coatings should consider the 
operating environments as a main design factor to maintain the tribological performance of 
coatings. Although a typical material has a limited operating range to maintain its best 
tribological performance, the nanocomposite coating design provides the multi-functionalities 
such as high hardness, high toughness, and low friction to one coating system for wide operating 
ranges. Previous developments on nanocomposite coatings show the improvement of hardness, 
toughness, friction reduction and wear resistance under severe operating conditions [1-8]. One 
of the nanocomposite design concepts is that nanocrystalline phases are encapsulated into an 
amorphous matrix because nanocrystalline phases with strong interaction with matrix phase 
provide superhardness. Veprek et al. [2] proposed a nanocomposite design concept based on 
avoiding dislocation fonnation and blocking the crack propagation in a 0.3 - 0.5 mn thin 
amorphous tissue. Also, Veprek and Argon [3] suggested an optimal nanocomposite grain size, 
3-5 mn, based on the interatomic bond distance. 

Yttria-stabilized-zirconia (YSZ) has a unique property combination that makes YSZ be 
an excellent candidate for a thin film to protect metal and alloy surfaces from high temperature 
environment, erosion and wear. However, since a ceramic is a very brittle material, YSZ films 
easily experience cracking when the deformation is beyond the elastic limit of the YSZ material. 
To improve the toughness and brittleness of oxide ceramic films, nanocomposite coatings with 
YSZ and gold were developed [4]. The combination of nanocrystalline phases and amorphous 
matrix improved the toughness while maintaining the high hardness and elastic modulus. Also, 
YSZ in a gold matrix was developed with encapsulated nanosized reservoirs of M 0 S 2 and 
diamondlike carbon (DLC) [1] because of the outstanding performance at high temperature. 
Furthermore, Voevodin et al. [1] proposed a “chameleon design concept” for 
YSZ/AU/M 0 S 2 /DLC and WC/DLC/WS nanocomposite coatings because of the unique coating 
skin adaptation, that is, “Chameleon” nanocomposite coatings show the reversible adaptability to 
humidity or temperature. “Chameleon” nanocomposite coatings were developed for the 
extremely changing environments from earth atmosphere to space by the self-adjusting surface 
tribological properties to maintain the perfonnance, that is, the coating friction surface changes 
its chemistry and structure depending on the operating environments. It was found that the 
mechanical properties of the nanocomposite coatings such as WC/WS 2 /DLC and 
YSZ/AU/M 0 S 2 /DLC exhibit a very strong dependence on coating composition [1]. 

Contact problems of the coated surfaces have been modeled by many researchers [9 - 14] 
using finite element analysis (FEA) technique. Most of the coatings investigated are single or 
multilayered coating systems. These studies were focused on the elastic or elastic-plastic 
behavior of the coating system and subsequent stress distribution under the Hertzian contact 
conditions. Parametric studies show the effects of critical design parameters such as contact 
width, friction coefficient, stiffness of layer, and etc. on the stress distribution and deformation 
inside the coating system. Coating parameter effects on the initiation and propagation of plastic 
zone for TiN/Ti coating system had been investigated using FEA models [9,10]. FEA model 
was developed to study the defonnation associated with indentation of coating system [11]. The 
elastic compliance of the coating system was detennined as a function of the indentation depth 
from the load-displacement curve. The stress-strain relationship between the substrate strength 


1 



and layer thickness was studied on TiN coating system and the consequent effect on the load 
capacity was investigated [12]. Gorishnyy and Olson [15] developed a FEA model to investigate 
the distribution of stress in single layer, bi-layer and multilayer films under combined normal 
and tangential loads. Recently, Veprek et al. [16] developed a new method for FEA modeling of 
ultra-hard ceramic coating indentations. Bull [17] repeated the FEA modeling of the scratch test 
for TiN on stainless steel with better substrate mechanical properties. 

This present study was focused on the study of stress distribution and elastic-plastic 
deformation on the nanocomposite coating systems developed at the Air Force Research 
Laboratory (RXBT). The FEA technique has been applied to the nanocomposite coating system 
model due to the nonlinear nature of the coating system behavior using commercial softwares 
such as ABAQUS/CAE and ABAQUS/Standard. The FEA model has been used to understand 
the mechanisms involved in tribological contacts for a wide range of material properties, 
compositions, and applied loads. Ultimately, this study will provide a design guideline for 
nanocomposite coatings. 

Nanocomposite Tribological Coating 

Voevodin et al. [1] illustrated the need for the nanocomposite coatings in extreme 
environments and developed the nanocomposite coating system shown in Figures 1 (a) and 1 (b). 
Figure 1(a) shows the schematic of nanocomposite coating system, featuring a 
crystalline/amorphous design in the main coating layer for cohesive toughness and a functionally 
gradient interface layer for adhesive toughness. It shows the nanocomposite coating, substrate, 
and gradient interface to enhance the adhesion strength and relieve interface stresses. Figure 1(b) 
depicts the schematic of nanocomposite coating design structure, that is, an amorphous matrix 
and crystalline phase (YSZ) were used for an optimum mechanical performance and load support. 
This nanocomposite has 3-10 mn crystalline grains in an amorphous matrix. The grains are 
separated by an amorphous matrix with 1-3 mn. This crystalline grain size is for the optimal 
hardness and toughness [1,3]. Also, this design was used to produce the optimum tribological 
performance. 

Figure 2 shows HRTEM (high resolution transmission electron microscope) image and 
SAD (submicron selected area electron diffraction) pattern obtained from (YSZ)o. 4 g Auo .12 
(M 0 S 2 V 18 C 0 . 22 , that is, an example of “Chameleon” nanocomposite coating produced by a 
hybrid of laser ablation and magnetron sputtering [1]. The SAD pattern of (YSZ ) 0 .48 Auo .12 
(M 0 S 2 V 18 C 0.22 nanocomposite coating demonstrates the rings for tetragonal ZrC >2 as well as face 
centered cubic Au and the diffused halo for amorphous carbon. Dark regions correspond to Au 
and YSZ nanocrystallines while bright regions correspond to amorphous matrix. The size of 
YSZ crystalline is 3-5 mn for the optimal load support and stability [3], The size of Au 
crystalline is 3 - 8 mn for the high temperature lubricant. The carbon and M 0 S 2 are amorphous 
for the ambient lubrication and the space lubrication, respectively. The initial model was 
proposed and explained by [5,6]. The dislocation formation is restricted in small grains (3 - 5nm) 
and the strength of nanocrystals approaches the ideal one that is of the order of about 10 % of the 
shear modulus. 

Finite Element Analysis Approach 

Figures 3 (a) and 3 (b) show the schematic of line contact for nanocomposite coating 
system and its corresponding finite element analysis model. The two-dimensional line contact 
finite element analysis (FEA) model is to analyze contact stress distribution in nanocomposite 
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coating system generated by a rigid cylindrical indenter. The commercial finite element 
softwares such as ABAQUS/CAE and ABAQUS/Standard were used to model the coating 
system under the plane-strain assumption. To reduce the computing time, the relatively coarse 
mesh was used for the region far away from the contact zone. Very fine mesh was used for the 
nanocomposite coating adjacent to the coating/rigid indenter interface near the contact zone in 
order to determine the stress distribution and deformation field accurately. Nearly 127671 nodes 
and 253331 elements were used to model the entire system. The three-node linear plain strain 
element (CPE3) was used. Coulomb friction was used to model the contact between the top 
coating surface and the rigid indenter. Due to the perfect bonding assumptions, the continuous 
displacement at the grain-matrix interfaces, coating interfaces, and coating/substrate interface is 
applied. Since the strain hardening is not considered, elastic perfectly plastic stress-strain 
relationships were used for the gradient coating materials and the substrate material. The onset 
of plastic yield begins by the von Mises yield criterion for the isotropic materials. Figure 3(b) 
demonstrates the von Mises stress distribution inside the nanocomposite coating system 
generated by a rigid indenter. Please note that the residual compressive stresses associated with 
the process of deposition are negligible because the nanocomposite coating thickness is a few 
micros. Nanocomposite coating system design incorporating functionally gradient interfaces 
with nanocomposites provide significant improvement of the tribological properties [7]. The 
design suggests the small size hard crystalline grains surrounded by an amorphous matrix to stop 
crack propagation. Grain boundary sliding improves the ductility of the nanocomposite coatings. 

Figures 4 (a) and 4 (b) show the two dimensional random microstructure geometries 
generated based on the concept of Voronoi tessellation. A Voronoi tessellation represents a cell 
structure constructed from a Poisson point process by introducing planar cell walls perpendicular 
to lines connecting neighboring points. This results in a set of convex polygons embedding the 
points and their domains of attraction, which completely fill up the underlying space. The 
concept of Voronoi tessellation has recently been extensively used in material science, especially 
to model random microstructures like aggregates of grains in polyscrystals, patterns of 
intergranular cracks, and composites [18]. Figures 4 (a) and 4 (b) have been generated by using 
MATLAB. The RAND command produces pseudo-random numbers [19]. The sequence of 
numbers generated is determined by the state of the generator. Figure 4 (a) was generated by 
resetting the generator to its 75 th state while figure 4(b) was generated by resetting the generator 
to its 50 th state. 

Figure 5 depicts the finite element analysis model developed for nanocomposite coating 
system nanoindentation. The coating system consists of nanocomposite coating, gradient 
interface, and substrate. Figure 5 clearly demonstrates the random geometrical structure of 
nanocomposite coating and functionally gradient interface generated using ABAQUS/CAE. The 
enlarged view shows the combination of crystalline phase/amorphous matrix structure. Table 1 
shows the mechanical material properties of YSZ/Au/MoSi/DLC composite coatings [1]. 

Coating hardness and modulus of elasticity were evaluated using nanoindentation. The results 
show that the nanocomposite composition and deposition temperature play important roles for 
improving the hardness and modulus of elasticity of the coating. The hardness of the YSZ/Au 
coatings decreases with the addition of M 0 S 2 and, to a less extent, carbon. 

Yield Strength Calculation Procedure 

To develop the elastic-perfectly plastic FEA model for nanocomposite coating system 
with functionally gradient Ti/TiC interface, the yield strengths of thin interface materials should 
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be obtained. Although there are no measured yield strength results from nanoindentation 
experiments, by the aid of the numerical nanoindentation simulation using a FEA model, the 
yield strengths of thin interfaces are able to be calculated. Therefore, an interactive procedure 
between FEA model simulations and nanoindentation experiments is required to calculate the 
yield strengths of the thin interface materials. It is assumed that the plastic deformation of an 
interface material is initiated when its von Mises stress generated inside the interface block 
becomes greater than yield strength of the interface material. 

The measured modulus of elasticity of the interface material is used as an input parameter 
to the FEA model for the nanoindentation simulation. Since yield strength of the interface 
material is required as an input parameter for the nanoindentation simulation, the initial yield 
strength of the interface material is estimated based on the measured hardness value of the 
material, i.e., one third of the hardness. Then, the numerical load vs. displacement curve is 
generated using the initial estimated yield strength and the measured modulus of elasticity. If the 
loading and unloading parts of the load vs. displacement curve show the discrepancy between the 
experimental and numerical curves, the second numerical load vs. displacement curve is needed 
to be obtained with either reduced or increased estimated yield strength. The slope of the second 
numerical load vs. displacement curve indicates the magnitude of the next estimated yield 
strength value. This step requires an iterative procedure to obtain a best-fit load vs. displacement 
curve from the nanoindentation simulation by changing the yield strength input parameter value 
for the simulation. Therefore, when the computational best-fit load vs. displacement curve 
matches the experimental averaged load vs. displacement curve closely, the estimated yield 
strength value is considered as the yield strength of the material. Please note that the yield 
strength of each material is determined from the nanoindentation simulation. 

Results and Discussion 

The experimental results for Ti .53C.47 interface from nano indentation are shown in 
Figures 6 (a) through 6 (d). The indentation depth to Ti.53C.47 interface sample is limited to 100 
mn to avoid the substrate effect on measured material properties. Figure 6(a) shows scattered 
data for hardness vs. displacement curves because of the surface roughness effect on the 
measurement. The hardness increases as the nanoindentation depth increases for Ti 53C.47 
interface. The average hardness increases as the nanoindentation depth increases for Ti.53C.47 
interface. The average hardness of Ti.53C.47 interface from the nanoindentation experiments is 
14.2 ± 1.8 GPa. The thick solid line represents the average hardness value of Ti.53C.47 interface 
along the indentation depth. Figure 6(b) depicts the modulus of elasticity vs. displacement 
curves for Ti.53C.47 interface. The average modulus of elasticity of Ti.53C.47 interface is 231.1 ± 
34.6 GPa. The thick solid line shows the average value of the modulus of elasticity for Ti.53C.47 
interface along the indentation depth. Figure 6(c) illustrates the load vs. displacement curves for 
Ti.53C.47 interface from nanoindentation. Due to the roughness of surface, the data is scattered. 
Figure 6(d) shows the comparison of averaged experimental results and best-fit FEA results for 
Ti.53C.47 interface. The average load vs. displacement curve agrees well with best-fit load vs. 
displacement curve when Y = 10 GPa and E = 240 GPa are used for Ti.53C.47 interface 
nanoindentation simulation. 

Figure 7 demonstrates the schematic of functionally gradient (FG) Ti/TiC interface 
design with programmed variation of compositional and mechanical properties proposed by 
Voevodin et al. [ 20 ]. Each layer has its unique material properties and thickness. Hardness and 
modulus of elasticity of Ti/TiC interface system are measured from interface samples. 25 
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indentation locations per sample were averaged to obtain the statistical values. The yield 
strength for each interface material is calculated from FEA nanoindentation simulation using the 
measured hardness and modulus of elasticity. These are the values used for inputs to FEA 
nanocomposite coating system model for this study. The material properties of each layer were 
linearly interpolated from the measured material properties data. To improve adhesion between 
the substrate and upper carbide layers, a-Ti is used as a bond layer. 

The FEA model was used to evaluate; (1) the effect of different crystalline phases such as 
YSZ and M 0 S 2 in amorphous YSZ/Au matrix on the von Mises stress distribution; (2) the effect 
of different amorphous YSZ/Au matrix compositions on the von Mises stress distribution as YSZ 
crystalline phase molar % increases; and (3) the effect of different crystalline phase densities 
such as 20 Vol. % and 75 Vol. % (Volume %) in amorphous YSZ/Au matrix on the von Mises 
stress distribution. Figures 8 (a) through 8 (d) show the contour plots of von Mises stress 
distribution inside the nanocomposite coating system indented by a rigid cylindrical indenter in 
order to study the effect of different crystalline phase properties on von Mises stress distribution 
inside the different nanocomposite coatings such as YSZ. 71 Au.nYSZ, YSZ /Au 11 M 0 S 2 , 
YSZ. 58 Au.i 8 YSZ, and YSZ/Au. isMoSt. Figures 8 (a) through 8 (d) illustrate that the von Mises 
stress field below the surface varies as the crystalline phase and YSZ/Au matrix composition are 
changed with indentation depth = 7 nm and friction coefficient, 0.0. Please note that the 
maximum von Mises stress always occurs at certain depth along the axis of depth below the 
center of the contact. Figure 8 (a) is the contour plot of von Mises stress for YSZ. 71 Au.nYSZ 
nanocomposite coating with a normal displacement, 8 = 7 nm. The maximum von Mises stress 
occurs inside the YSZ. 71 Au.nYSZ coating system with the magnitude of 2.28 GPa. The YSZ 
crystalline phase exhibits the highest von Mises stress because of the higher modulus of elasticity 
than YSZ. 7 iAu.ii matrix. With the M 0 S 2 crystalline phase, the maximum von Mises stress occurs 
below the contact surface with the magnitude of 1.90 GPa inside the YSZ. 71 Au .11 amorphous 
matrix. Figure 8 (c) shows the contour plot of von Mises stress for YSZ. 58 Au.i 8 YSZ with a 
nonnal displacement, 8 = 7 nm. The maximum von Mises stress occurs with the magnitude of 
2.27 GPa. Figure 8 (d) shows the contour plot of von Mises stress for YSZ. 58 Au.i 8 MoS 2 with a 
nonnal displacement, 8 = 7 nm. The maximum von Mises stress occurs with the magnitude of 
1.90 GPa. Note that there are vertical lines on the contour plot of von Mises stress distribution 
inside the nanocomposite coating because of the repeated duplication of the small nanocomposite 
structure pattern. In general, YSZ. 71 Au.nYSZ nanocomposite coating generates higher von 
Mises stress distribution than YSZ 71 AU. 11 M 0 S 2 nanocomposite coating under the same loading 
conditions because of the higher Young’s modulus of YSZ crystalline phase. YSZ. 71 Au.nYSZ 
nanocomposite coating has slightly higher von Mises stress than YSZ 58 Au.i 8 YSZ nanocomposite 
coating because YSZ. 71 Au .11 amorphous matrix has higher Young’s modulus than YSZ. 58 Au.i 8 
matrix. Table 2 shows the material property values used for the FEA model [ 8 ]. 

Figures 9 (a) through 9 (d) show the contour plots of the von Mises stress distribution for 
YSZ/AU/M 0 S 2 nanocomposite coatings with different amorphous matrix compositions to 
represent the effect of amorphous matrix composition properties on Von Mises stress distribution. 
Figure 9(a) depicts the contour plot of von Mises stress distribution inside YSZ. 65 Au. 35 M 0 S 2 
nanocomposite coating system with the indentation depth of 7 nm with maximum von Mises 
Omax = 1-27 GPa. The maximum von Mises stress occurs YSZ/Au amorphous matrix region 
because of the higher Young’s modulus of YSZ/Au matrix. Figures 9(b) through 9(d) show the 
contour plot of von Mises stress distribution inside YSZ. 75 AU. 25 M 0 S 2 nanocomposite coating 
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with the maximum von Mises stress, a max = 1.83 GPa, YSZ/Au. 15 M 0 S 2 with a max = 2.22 GPa, 
and YSZ. 95 AU. 5 M 0 S 2 with a max = 2.59 GPa, respectively. In general, the hardness and young’s 
modulus of YSZ/Au matrix increases as the YSZ composition increases in YSZ/Au matrix. 
Therefore, the YSZ. 95 AU. 5 M 0 S 2 nanocomposite coating generated the highest maximum von 
Mises stress when the displacement of the indenter tip, 7 nm, is applied to all the cases. The 
table 3 shows the material properties used for the FEA model [8]. 

Figures 10 (a) through 10 (b) show the contour plots of equivalent plastic strain for 
different nanocomposite coatings after the nanoindenter tip is removed from the nanocomposite 
coating surface. Figure 10 (a) illustrates contour plot of the equivalent plastic strain distribution 
generated inside the YSZ/Au. 35 M 0 S 2 nanocomposite coating. It indicates that the maximum 
equivalent plain strain occurs in the M 0 S 2 with 8.46 x 10" 4 . Figure 10 (b) depicts the contour 
plot of equivalent plastic strain distribution generated inside the YSZ. 75 AU. 25 M 0 S 2 
nanocomposite coating. The magnitude of the equivalent plastic strain increases up to 2.30 x 10" 
3 in M 0 S 2 nanocrystalline phase. The area of the pennanent deformation increases because the 
higher modulus of elasticity of YSZ/Au matrix causes higher von Mises stress in 
YSZ. 75 AU. 25 M 0 S 2 nanocomposite coating. Figure 10 (c) demonstrates the contour plot of 
equivalent plastic strain distribution generated inside the YSZ/Au. 15 M 0 S 2 nanocomposite 
coating. The magnitude of the equivalent plastic strain increases up to 2.91 x 10" . Figure 10 (d) 
displays the contour plot of equivalent plastic strain distribution generated for YSZ/Au. 5 MoS 2 
nanocomposite coating. In general, as the modulus of elasticity increases, the magnitude of the 
von Mises stress generated increases. Therefore, the magnitude of the equivalent plastic strain 

a 

increases up to 3.48 x 10" .because of the higher modulus of elasticity of YSZ/Au. 5 M 0 S 2 
nanocomposite coating. In addition, as the YSZ composition increase in YSZ/Au amorphous 
matrix, the larger area below the contact is affected for the plastic deformation, that is, larger 
pennanent plastic deformation area occurs because of the modulus of elasticity and hardness 
increase. 

During the nano indentation process, the YSZ/AU/M 0 S 2 or YSZ/Au/YSZ coating system 
still experience the residual stress after the nanoindenter tip is removed from the nanocomposite 
coating system. Figures 11(a) through 11(d) show the contour plot of residual von Mises stress 
distribution generated inside the nanocomposite system obtained from the FEA model. As the 
YSZ composition increases in YSZ/Au matrix, the nanocomposite coating system shows the 
higher residual von Mises stresses generated inside the coating. Figure 11 (a) depicts the contour 
plot of residual von Mises stress in YSZ/Au. 35 M 0 S 2 nanocomposite coating after the 
nanoindentation depth of 7 nm applied and removed from the coating surface. The maximum 
von Mises stress, 7.26 X 10 7 Pa, is generated below the center of the contact. Figure 11 (b) 
demonstrates the contour plot of residual von Mises stress distribution inside the 
YSZ. 75 AU. 25 M 0 S 2 nanocomposite coating. The maximum residual von Mises stress generated 
below the contact is 2.22 X 10 Pa. As the YSZ composition increases in YSZ/Au amorphous 
matrix, the larger region below the contact surface is influenced by the residual stress under the 
same loading conditions. Figure 11 (c) illustrates the contour plot of residual stress distribution 
inside the YSZ/Au. 15 M 0 S 2 nanocomposite coating. The maximum von Mises stress generated 
below the contact is 3.22 X 10 Pa. Figure 11 (d) displays the contour plot of residual von Mises 
stress distribution inside the YSZ. 95 AU. 5 MoS 2 nanocomposite coating. The maximum residual 
von Mises stress generated below the contact is 4.17 X 10 8 Pa. Larger area of the coating system 
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has the residual von Mises stress inside the nanocomposite coating as the YSZ composition 
increases. 

Figures 12 (a) through 12 (d) show the contour plots of von Mises stress distribution 
inside the nanocomposite coating system to represent the effect of density of crystalline phase 
material composition on Von Mises stress. Figure 12 (a) depicts the von Mises stress generated 
inside YSZ. 75 AU 25 M 0 S 2 nanocomposite coating with 75 phase volume % density. Figure 12 (b) 
displays the von Mises stress generated inside YSZ. 75 AU. 25 M 0 S 2 nanocomposite coating with 20 
phase volume % density. Because of the lower density (20 %) of the M 0 S 2 crystalline phase, the 
higher maximum von Mises stress, a max = 3.67 GPa, was generated. Figure 12 (c) demonstrates 
the von Mises stress distribution generated inside the YSZ. 95 AU. 5 M 0 S 2 nanocomposite with 75 
phase volume % density. Figure 12 (d) illustrates the von Mises stress distribution generated 
inside the YSZ. 95 AU. 5 M 0 S 2 nanocomposite with 20 phase volume % density. Again, higher 
maximum von Mises stress, a max = 5.43 GPa, was generated in lower density (20 %) 
nanocomposite coating. In general, higher density (75 %) nanocomposite coating generates 
lower maximum von Mises stress while the higher density nanocomposite can generate a narrow 
gap between crystalline grains. Therefore higher volume density nanocomposite is desirable for 
the same loading conditions. 

Figures 13 (a) and 13 (b) show the contour plots of von Mises stress distribution to 
represent the effect of indenter size on Von Mises stress inside YSZ. 75 AU. 25 M 0 S 2 nanocomposite 
coating system. Figure 13 (a) depicts the von Mises stress distribution in YSZ. 75 AU. 25 M 0 S 2 
nanocomposite coating system with indenter radius, 7.5 pm, and indenter displacement, 6 nm. 
The maximum von Mises stress generated is 1.71 GPa under the contact surface. Figure 13 (b) 
illustrates the von Mises stress in YSZ. 75 AU. 25 M 0 S 2 nanocomposite coating system with indenter 
radius, 15 pm, and displacement, 8 nm. The maximum von Mises stress occurs is 1.64 GPa. 
Since the area of von Mises stress generated by indentation is significantly influenced by the 
indenter size, the nanocomposite coating thickness design requires a careful consideration on the 
indenter size. 

Figures 14 (a) through 14 (b) show the contour plots of von Mises stress distribution 
inside the YSZjjAu. 11M0S2 nanocomposite coating to represent the effect of applied load on von 
Mises stress. The different indentation depths from 5 nm to 11 nm applied to the 
nanoindentation simulation for the same contact conditions. Figures 14 (a) and 14 (b) depict the 
von Mises stress distributions inside the YSZjiAu. 11M0S2 nanocomposite coating for indentation 
depth, 8 = 5 nm and 7 nm. Higher von Mises stress is generated below the contact surface due to 
the higher displacement, 7nm. Figure 14 (c) shows the von Mises stress inside the 
YSZ 71 Au 11M0S2 nanocomposite coating for indentation depth, 8 = 9 nm. Maximum von Mises 
stress increases up to 2.96 GPa due to the increased indentation depth. Figure 14 (d) depicts the 
von Mises stress distribution of YSZ71AU11M0S2 nanocomposite coating with 8=11 nm. In 
general, as the indenter displacement increases, wider Hertzian contact width and higher von 
Mises stress occur. 

Summary and Conclusion 

A FEA model has been developed to investigate the stress distribution and elastic-plastic 
deformation generated inside the nanocomposite coating system due to the nanoindentation. The 
concept of Voronoi tessellation was applied to generate the random microgeometries of 
nanocomposite coatings. In general, YSZ/Au amorphous matrix provided excellent mechanical 
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properties for contact load support. YSZ/Au coatings with the addition of YSZ crystalline have 
better load supporting mechanism than YSZ/Au coatings with the addition of M 0 S 2 . The 
increase of the YSZ composition in YSZ/Au amorphous matrix caused an increase in elastic 
modulus and subsequent load supporting. Also, the increase of the YSZ composition in YSZ/Au 
amorphous matrix caused an increase in equivalent plastic strain and subsequent residual stress 
under the same loading conditions. However, the magnitude of the residual von Mises stress is 
relatively smaller than the stress generated during the deposition process. The volume density of 
the crystalline phase material composition plays an important role to maintain the integrity of the 
coating system because of the material properties and the gap between the crystalline phases. 
Also, the indenter size effect is significant for nanocomposite coating thickness design. In 
addition, the applied load is the important operating parameter to design the nanocomposite 
coating system because the applied load determines coating thickness, crystalline phase density, 
YSZ/Au composition, and etc. 
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Table Captions 

Table 1. Mechanical material properties of YSZ/AU/M 0 S 2 /DLC composite coatings 
produced at different deposition temperature [1]. 

Table 2. Mechanical Material properties of YSZ/Au matrix, YSZ, and MoS2. 

Table 3. Mechanical Material properties of YSZ/Au matrix with different amorphous 
matrix composition properties. 
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Table 1. Mechanical material properties of YSZ/AU/M 0 S 2 /DLC composite coatings 
produced at different deposition temperature [1], 


Coating Designation using molar 
fractions of the components 

Deposition 
temperature (°C) 

Modulus (GPa) 

Hardness (GPa) 

(YSZ)o.82(MoS2)o.18 

150 

118 

10 

(YSZ)o.7iAuo.ii(MoS2)o.i8 

150 

47 

2 

(YSZ)o.58Auo.18(MoS2)o.24 

300 

131 

4 

(YSZ)o.9lCo.09 

150 

173 

15 

(Y SZ)o.78AUo. 13C0.09 

150 

66 

11 

(Y SZ) 0 .8 1 Auo. 15C0.04 

300 

187 

13 

(YSZ)o.45(MoS2)o.33C 0 .22 

150 

56 

5 

(YSZ)o.62(MoS2)o.2lCo.l7 

300 

169 

12 


Table 2. Material properties of YSZ/Au matrix, YSZ, and MoS2. 


Material 

Modulus of Elasticity [Pa] 

Yield Strength [Pa] 

YSZ. 7 iAu.ii 

2.4 x 10 11 

5.5 x 10 9 

YSZ.58Au.l8 

2.0 x 10 11 

3.7 x 10 9 

YSZ 

3.2 x 10 11 

8.3 x 10 9 

M 0 S 2 

1.2 x 10 11 

1.0 x 10 9 


Table 3. Material properties of YSZ/Au matrix with different amorphous matrix 

composition properties. 


Material 

Modulus of Elasticity [Pa] 

Yield Strength [Pa] 

YSZ. 95 AU .5 

2.6 x 10 11 

7.3 x 10 9 

YSZ. 85 AU .15 

2.3 x 10 11 

5.3 x 10 9 

YSZ. 75 Au .25 

2.0 x 10 11 

3.7 x 10 9 

YSZ. 65 Au 3 5 

1.5 x 10 11 

2.8 x 10 9 
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adaptive transfer film (“tribo-skin”) on 
contact surfaces 



Substrate 


wear debris 

“chameleon”coating 
with lubricant reservoirs 


gradient interface 


(a) 



3-10 nm 


1-3 nm 


solid lubricant 
nanoparticle 


\ hard crystalline 

nanoparticle 


(b) 


Figure 1. (a) Schematic of nanocomposite coating system, featuring a 
crystalline/amorphous design in the main coating layer for cohesive toughness and a 
functionally gradient interface layer for adhesive toughness; (b) Schematic of 
nanocomposite coating design structure, that is, an amorphous matrix and crystalline phase 
(YSZ) were used for an optimum mechanical performance and load support. 
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Figure 2 . HRTEM (high resolution transmission electron microscope) image and 
SAD (submicron selected area electron diffraction) pattern obtained from (YSZ)o.48 
Auo.12 (MoS2)o.i8 C0.22, that is, an example of “Chameleon” nanostructure, i.e. 
YSZ/AU/M0S2/DLC nanocomposite coating produced by a hybrid of laser ablation 
and magnetron sputtering [ 1 ], 




w 


Line Contact 



Nanocomposite 

Coating 

Gradient 

Interface 

Substrate 



(a) 


(b) 


Figure 3. (a) Schematic of line contact for nanocomposite coating system; (b) Finite 
element analysis model for nanoindentation of nanocomposite coating system. 
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(a) 



(b) 


Figure 4. Schematic of random microstructure geometries generated 
based on the concept of Voronoi tessellation; (a) 75 th state; (b) 50 th 
state. 
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Combination of nanocrystalline/ 


amorphous structure 
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Figure 5. Finite element analysis model of nanocomposite coating system 
nanoindentation. 
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Figure 6. Nanoindentation results for Ti.53C.47 coating; (a) hardness vs. displacement plot, (b) modulus of elasticity vs 

displacement plot, (c) load vs. displacement plot during loading and unloading stages, and (d) comparison of measured 
average load vs. displacement curve and best-fit load vs. displacement curve from experiments and FEA simulations, 
respectively. 
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Figure 7 Schematic of functionally gradient Ti/TiC Interface Design 


19 



















































oi cri ot Q> Qi o» at co oo oo oo oo oo d 
oooooooooooooq 
' + + + + + + + + + + + + + • 

' OJ ILI OJ <U Cl <U OJ OJ >U <U <U 0.1 f 

. t"-o cm co m ["-oo o r-co o ["-co c 
r'O'rcorg'ooortcoin'ccot 

cm cti r-in tr <m»h in cri co r-^-iinfc 




ca o o 
£ U U 




ca 

o 

U 



20 


Figure 8. Contour plots of von Mises stresses distribution generated inside the nanocomposite coating system from the FEA 
model with different crystalline properties at indentation depth, 8 = 7 nm: (a) YSZ.71Au.nYSZ; (b) YSZ.71AU.11M0S2; (c) 
YSZ.jgAu.igYSZ; and (d) YSZ.sgAu.igMoSi- 
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Figure 9. Contour plots of von Mises stresses distribution generated inside the nanocomposite coating system from the FEA 
model with different amorphous matrix composition properties at indentation depth, 8 = 7 nm: (a) YSZ. 65 Au 3 5 MoS 2 , a max 
=1.27 GPa; (b) YSZ. 75 Au. 25 MoS 2 , a max = 1.83 GPa; (c) YSZ 85 AU. i 5 MoS 2 , a max = 2.22 GPa; and (d) YSZ. 95 Au 5 MoS 2 , a max = 
2.59 GPa. 
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Figure 10. Contour plots of equivalent plastic strain distribution generated inside the nanocomposite coating system from the 
FEA model with different amorphous matrix composition properties at indentation depth, 8 = 7 mn: (a) YSZ 65 Au. 35 MoS 2 , 8.46 
x 10' 4 ; (b) YSZ. 75 AU. 25 M 0 S 2 , 2.30 x 10' 3 ; (c) YSZ. 85 AU. 15 M 0 S 2 , 2.91 x 10' 3 ; (d) YSZ 95 Au 5 MoS 2 , 3.48 x 10' 3 . 
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Figure 11. Contour plots of residual von Mises stress distribution generated inside the nanocomposite coating system from the 
FEA model with different amorphous matrix composition properties at indentation depth, 8 = 7 nm: (a) YSZ 65 Au. 35 MoS 2 , 7.26 
x 10 7 Pa; (b) YSZ. 75 AU. 25 M 0 S 2 , 2.216 x 10 8 Pa ; (c) YSZ 85 AU. 15M0S2, 3.22 x 10 8 Pa; (d) YSZ. 95 Au 5 MoS 2 , 4.17 x 10 8 Pa. 
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Figure 12. Contour plots of von Mises stress distribution generated inside the nanocomposite coating system from the FEA 
model with different matrix/inclusion material composition at indentation depth, 8 = 7 mn: (a) YSZ.7 5 Au.25MoS2, a max = l-83 
GPa; (b) YSZ.75AU.25M0S2, a max = 3.67 GPa; (c) YSZ. 95AU5M0S2, a max = 2.59 GPa; (d) YSZ.95AU.5M0S2, a max = 5.43 GPa. 
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Figure 13. Contour plots of von Mises stress distribution generated inside the 
nanocomposite coating system from the FEA model with different indenter radius and 
indentation depth (a) YSZ.75AU.25M0S2, R = 7.5 pm, 8 = 6 nm; (b) YSZ.75AU.25M0S2, R = 
15 pm, 8 = 8 nm. 
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Figure 14. Contour plots of von Mises stress distribution generated inside the YSZ 7 iAu 11 M 0 S 2 nanocomposite coating system 
from the FEA model at different indentation depths, 8 = 7 mn: (a) 5 = 5 nm; (b) 8 = 7 nm; (c) 8 = 9 nm; (d) 8 = 11 mn. 































